Models of quintessence coupled to the electromagnetic field and the cosmological 

evolution of alpha 
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We study the change of the effective fine structure constant in the cosmological models of a scalar 
field with a non-vanishing coupling to the electromagnetic field. Combining cosmological data and 
terrestrial observations we place empirical constraints on the size of the possible coupling and explore 
a large class of models that exhibit tracking behavior. The change of the fine structure constant 
implied by the quasar absorption spectra together with the requirement of tracking behavior impose 
a lower bound of the size of this coupling. Furthermore, the transition to the quintessence regime 
implies a narrow window for this coupling around lO^'' in units of the inverse Planck mass. We 
also propose a non-minimal coupling between electromagnetism and quintessence which has the 
effect of leading only to changes of alpha determined from atomic physics phenomena, but leaving 
no observable consequences through nuclear physics effects. In doing so we are able to reconcile 
the claimed cosmological evidence for a changing fine structure constant with the tight constraints 
emerging from the Oklo natural nuclear reactor. 
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INTRODUCTION 



Recently claimed observational evidence for a ^ lOppm 
change of the fine structure constant over cosmological 
time between z ~ 0.5 and 2 = brings an unexpected 
and interesting support to the old idea of Dirac Jj . The 
result of Webb et al. |2| (if it continues to defy alternative 
explanations by systematic effects ^) has serious impli- 
cations for cosmology, astrophysics and particle physics, 

The most powerful conclusion that follows from [3| is 
the existence of new fields, massless or nearly massless, 
and different from gravitation, electromagnetism, neutri- 
nos and axions. Indeed, the "promotion" of a coupling 
constant g into a function of time immediately implies 
that it should also depend on the rest of the coordinates, 
g = g{x). The dependence of g on a: in an interacting 
theory will inevitably generate a kinetic term d^gd^g, 
so that the "coupling constant" becomes a propagating 
field. In order to avoid confusing terminology of "chang- 
ing alpha" , it is reasonable to introduce a new degree of 
freedom, a scalar field (j){x), via a chain of relations: 



J_ TP pfil' 



— BF{(b{x))F^,F^' 
% 



(1.1) 



We postulate that the scalar field couples to the electro- 
magnetic F^^F^^'^ via some unknown dimensionlcss func- 
tion Bp((j){x)). Using the freedom in the definition, one 
can always choose Bp{(j){z = 0)) = 1, so that go is the 
coupling constant "today" , i.e. at the redshift z = 0. 



The dependence of gauge coupling constants on coor- 
dinates can hardly be regarded as a theoretical novelty. 
Indeed, in string theory/supergravity all effective cou- 
pling constants depend on scalar field(s) (j)i{x), so that 
g — g{4)i{x)). However, the late change of g, a few 
billion years after the Big Bang, implies that some of 
these fields have essentially near fiat potentials, i.e. their 
masses must be comparable to the Hubble parameter at 
the present epoch ~ 10~^^ eV. This is in contrast with 
the phenomenologically developed models of string the- 
ory/supergravity, where all moduli are believed to receive 
electroweak scale masses after supersymmetry breaking, 
which makes them completely frozen long before nucle- 
osynthesis, not mentioning the present epoch. Therefore, 
such a late evolution of the effective a comes as a total 
surprise for particle physics. 

In general, a nearly massless scalar poses a serious 
problem of naturalness, which only worsens if this scalar 
is interacting fl^. Indeed, the inclusion of a non- 
renormalizable interaction of the form Bp{<j))F'^ at the 
quantum level requires the existence of a momentum cut- 
off A j/y. Unfortunately,, any particle physics motivated 
choice of Ajjv destabilizes the quintessence potential, i.e. 
it could induce a mass term much larger than the re- 
quired ^ Hq. However, because the nature of this 
fine-tuning is similar to the one required for the small- 
ncss of the cosmological constant, we choose to proceed 
hoping that one day both problems could be resolved si- 
multaneously. 

The study of the interaction has a long history. 
Initially, it was examined by Bckenstein 4], who intro- 
duced the exponential form for the coupling of the scalar 
field to the electromagnetic Lagrangian which in prac- 
tice can always be taken in the linear form Bp{(j}{x)) — 
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1 — C'^0> where = 8ttG. He showed that the exchange 
by this scalar leads to the effective non-universality of 
the gravitational force which bounds the possible size of 
C at one per mil level. The cosmological evolution of 
is dynamical and it starts immediately after the transi- 
tion to a matter-dominated Universe. In principle, the 
existence of a non-zero ^ alone is sufficient to ensure the 
cosmological evolution of 0, driven by the electromag- 
netic portion of the baryon mass density. However, in 
the minimal Bekenstein model where ^ is the only source 
for (/), the resulting change of a is way too small to be 
interesting: \a{z - 1) - ao|/ao ^ 10"^° - lO'^ HHIIl 
An enhancement by up to five-six orders of magnitude 
may occur by abandoning the minimalistic approach and 
introducing a "generalized" Bekenstein model where 
is driven either by 0(1) couplings to dark matter or by 
its own potential 0, 0|, in which case Aa/a may be 
consistent with the results of Refs. 

Quite independently from the physics of changing cou- 
plings, there has been a tremendous increase of interest 
towards cosmological models with a light scalar field as 
a possible source of dark energy or quintessence. Evi- 
dence for the change of the effective coupling constant is 
a strong argument in favor of the existence of a scalar 
field with the wave length comparable to the size of 
the Universe which thus makes it a good candidate for 
quintessence. Even though the concept of quintessence 
and the choice if its potential looks completely ad hoc, 
some quintessence models have an advantage over a "sim- 
ple" explanation of dark energy by a non-zero cosmologi- 
cal constant. Indeed, certain forms of potentials allow for 
an attractor-type solution, in which case the late cosmo- 
logical evolution of has no memory of the initial con- 
ditions that could be chosen (almost) arbitrarily. The 
evolution of the energy density of the (fi field is subdom- 
inant with respect to the energy density of matter and 
radiation, and only during the last stages of the cosmo- 
logical evolution does it become a dominant component. 
It is thus intriguing to identify the field (p responsible 
for the time dependence of a with tracking quintessence 
and study possible patterns for Aa/a as a function of 
cosmological time. 

In this paper, we present a thorough analysis of a large 
class of scalar field models with tracking behavior sup- 
plemented by a coupling to the electromagnetic field. A 
similar, but more restricted comparison has recently been 
made in [s^ ■ Where appropriate we will compare our re- 
sults to those. 

A successful model must satisfy the following present 
constraints: 



(1) Oklo natural reactor (z = 0.14) 
Aa 



< 10" 



-7 . 



(2) Decay rate ^^^Re ^^^Os, (z = 0.45) jH: 

<10-; 



(1.2) 



(1.3) 



(3) Cosmic microwave background radiation (z = 10"^) 
133: 



Ac 



< 10" 



(1.4) 



(4) Big Bang nucleosynthesis (z = 10*^ - 10^°) [il 
<2xl0-^ 



1.5) 



(5) Atomic clocks (z = 0) 

< 4.2 X 10"^^ yr~^ , 



(1.6) 



where a dot represents differentiation with respect 
to cosmic time; 



(6) Equivalence principle tests [Tj|: 

Id < 10-3. 



(1.7) 



In addition, we take into account the claimed evidence 
for the variation of the fine structure constant. We will 
assume 



(7) Quasar absorption lines (z = 3) 0|: 
Aa 



= -10 



-5 



(1.8) 



We first address a minimal Bekenstein-like coupling 
and find the behavior of the fine structure constant as 
a function of the redshift for the exponential, power-like 
and exponential times power-like potentials. 

In all cases we analyze the scalar field evolution in the 
tracking regime i.e., when the energy density of the scalar 
field remains sub dominant. In this case we are able to 
present the lower bound on consistent with tracking 
behavior and a non-zero change of the fine structure con- 
stant The upper limit in this case is furnished by 
the direct experimental limit on the non-universality of 
the gravitational force. As a separate exercise, we con- 
sider a scalar field as quintessence, i.e. when the tracking 
regime changes to the scalar field domination at z ~ 1. 
In this case, we are able to determine the allowed win- 
dow for ^, without referring to the constraints from the 
non-universal gravitational forces. Also in the context of 
a quintessence model, we consider the case of the tracker 
field driven by the combination of the self potential and 
the coupling to dark matter. 

We confront the predictions for Aa/a with the above 
limits from various astronomical and terrestrial observa- 
tions, and find that the present day limits on a/a and 
constraints from the Big Bang nucleosynthesis are gen- 
erally safe. On the other hand, the Aa/a coming from 
the Oklo natural reactor [3, IS S I^S and from the mete- 
oritic abundance data |25j | are generally inconsistent with 



3 



tracking behavior. In the case of a non-minimal coupling 
of quintessence to the electromagnetic field, the results 
lose their predictivity as almost any profile for a{t) can 
be achieved by an appropriate choice of Bf{4>)- 

We go on to discuss a novel possibility that the effec- 
tive change of the coupling constant occurs only in atomic 
interactions while Aa/a remains effectively zero for nu- 
clear physics phenomena. The construction we propose 
is based on the idea that the coupling between photons 
and the scalar field may change for off-shell photons and 
has a form factor at some energy above the electron mass. 
This leads to the suppression of the coupling between the 
scalar field and the nuclear levels giving a chance to rec- 
oncile the results of Webb et al. and the bounds coming 
from Oklo phenomenon without fine-tuning. 

This paper is organized as follows: in the next sec- 
tion we analyze in detail four types of tracking scalar 
fields and determine the allowed values for C- Section 3 
addresses the issue of the non-minimal Bf{4>) and intro- 
duces the model where C has a form factor. We reach our 
conclusions is section 4. 



2. MINIMAL COUPLING OF QUINTESSENCE 
TO THE ELECTROMAGNETIC FIELD 

We study a class of models of a neutral scalar field 
coupled to electromagnetism: 



S= - 



1 

2^ 



d'^x\/~g R 

[C^ + Cb + C^f) , (2.9) 



where R is the scalar curvature. In this paper, we 
deal with the cosmological evolution and take the stress- 
energy tensor of the matter fields in a homogeneous form 
corresponding to a background perfect fluid with the 
energy-momentum tensor given by 



= {pB+PBy^-u" -g^.pB- (2.10) 

Throughout the paper, we employ (-1-,—,—,—,) signa- 
ture for the metric tensor. The pressure and en- 
ergy density are related by the equation of state 
— {j~l)pg. 7 is taken to be a constant (e.g. 7 — 4/3 
for radiation and 7 = 1 for matter) . The Lagrangian den- 
sity for the scalar field cf) is 



(2.11) 



The associated energy-momentum tensor of the scalar 
field then reads 



T^^, = d^'Pd.^P - g^, - y(0) 



(2.12) 



The interaction term between the scalar field and the 
electromagnetic field is 



(2.13) 



where Bf{4>) allows for the evolution in (/). In this 
paper we adopt a linear dependence on (f) such that 
Bf{4>) = 1 — C,i^{4> — 4'o), where the subscript repre- 
sents the present value of the quantity. The effective 
fine structure constant depends on the value of (j) as 
a = a^l Bf{4i)- Therefore, we have 



Ac 



a — ao 



Ck(0-(?!)o)- 



(2.14) 



We will consider a spatially flat Friedmann-Robertson- 
Walker (FRW) Universe with metric ds^ — dt^~o?{t)d:x? . 
We further assume a nearly homogeneous scalar field so 
that the spatial derivatives of the field can be neglected. 
Hence, we can define analogously to the case of a perfect 
fluid, the energy density and pressure of a scalar field as, 
P(t> = </'^/2 + V and p^ = </i^/2 — V, respectively. 
The governing equations of motion are 



4( 



IPb 



H = 

PB = -^iHpg, 



subject to the Friedmann constraint 



H' = ^iPB+P,)- 



(2.15) 
(2.16) 
(2.17) 

(2.18) 



The term C^f can be neglected in the equations deter- 
mining the dynamics of as its average is zero f or p ho- 
tons, and also quite small, O(10~'^) for baryons A 
notable exception is when couples strongly and directly 
to dark matter [l3|. 

In the following subsections we will consider different 
forms of the scalar potential V{(j)) and analyze two pos- 
sible cases of evolution of the field: (i) slow roll down the 
potential V{(f)) controlled by the friction provided by the 
background fiuid (tracker evolution); (ii) the field pro- 
vides the source of dark energy, and is just starting to 
dominate the dynamics leading the universe into a pe- 
riod of accelerated expansion. We will start by estimat- 
ing the value of the coupling ( consistent with Eq. H1.8|) 
by setting 



c 



10- 



3) - K(/)(0 = 0) 



(2.19) 



One should bare in mind that we are not attempting to 
fit the data to specific scalar potentials but instead to 
extract an estimate of C consistent with the data. 

We calculate the ratio of the variation of a with time 
at present using the relation 



-c 



d{K<p) 
d{l + z\ 



Hn 



(2.20) 
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where Hq is the value of the Hubble constant today, Hq = 
h/9.78 X 10"^ yr-i. We will use h = 0.65 throughout 
the paper. 



A. Exponential potentials 

For a pure exponential potential V{(t)) — 
Vq exp(— Ak(/)), there exist just two possible late 
time attractor solutions with quite different properties, 
depending on the values of A and the background's 
equation of state 7 [34|: 

(1) A^ > 37. The late time attractor is one where the 
scalar field mimics the evolution of the barotropic fluid 
with 70 = 7, and the relation fl^ = 87/ A^ holds. 

(2) A^ < 87. The late time attractor is the scalar field 
dominated solution (fJ^ = 1) with 7^ = A^/8. 

The second regime can be discarded as unrealistic. In- 
deed, there is a tight constraint on the allowed mag- 
nitude of at the time of Big Bang nucleosynthesis, 
MeV) < 0.045 Moreover, we must allow time 
for the formation of structure before the universe starts 
accelerating, which requires fl^ to be subdominant after 
the CMB decoupling. Therefore, the transition to fi^ ~ 1 
should occur at the present epoch. For this model, this 
is possible only at the expense of extreme fine tuning 
of the initial value p™, a type of fine-tuning that goes 
against the whole idea of the tracker/quintessence field. 
The first regime is perfectly compatible with all experi- 
mental data, if there is an additional component to dark 
energy which ensures the late time acceleration, as the (f> 
field itself in this model does not allow for an accelerating 
expansion. 

For now we simply assume that the scalar field respon- 
sible for the change in the value of a is tracking and is 
subdominant today, i.e. Jl^ < 1/2. This imposes a lower 
bound |A| > ^/E. The evolution of the field in the tracker 
regime is well known and is given by js^ 



-^ln(l + z) 
A 



(2.21) 



Therefore, combining Eqs. (|2.19|) and 12.21|1 . we obtain 



c 



10 



-5 



ln(4) 87 



(2.22) 



As we discussed before, if is subdominant today, then 
there must be some other source of dark energy respon- 
sible for the accelerating universe. The presence of a cos- 
mological constant introduces an error of less than 80% 
in the above estimate. With |A| being limited from be- 
low, this imposes the lower bound on the possible size of 
the coupling, C > 6 x 10"^. On the other hand, C is lim- 
ited from tests of the equivalence principle, |C| < 10~^, 
which we can use to further constrain the value of the 
parameter A in a matter dominated universe: 



V6<\X\< 415. 



From Eq. (|2.20|l we immediately extract that for this 
potential a/a = 4.8 x 10^^^ yr~^ consistent with the 
atomic clocks bound Eq. (|1.6(l . 

It is also instructive to check whether this range of 
couplings is consistent with constraints on Aa/a com- 
ing from the Big Bang nucleosynthesis. In the radiation 
dominated era, the evolution of the scalar field is given 

by 



K{(j) - 0o) 



8 4 
- ln(l -I- Zoq) - 



1 



, (2.24) 



where 1 -I- Zcq is the redshift at which the energy den- 
sity of the radiation and matter components have the 
same value, Zeq = ^m/^r ~ 5 x 10^. Substituting into 
Eq. (|2.14(l for the estimated value of (, it is easy to verify 
that Aa/a = -2 x 10"^ at redshift z = 10^°, which is 
fully consistent with the bound in Eq. H1.5|l . 



B. Power law potentials 

For an inverse power law potential of the form V{(f>) = 
M^+"/0" the cosmological evolution of the field is given 
by 36] 



- 0o) = K(t>o \-l + (1 + z)-'^^/(2+")' 



(2.25) 



where the value of at z = is given by 



8 K 1 jn 



7n 



-l/(2+n) 



2M4+"n2n + 2 V n + 

(2.26) 

Combining Eq. (|2.25|l with Eq. H2.19|l we can estimate 
the coupling with electromagnetism through 



c 



10- 



-1 + 4/3 



(k0o) ^ 



(2.27) 



(2.23) 



where /3 — —87/(2 + n). 

In Fig. n] we show a contour plot for log ^ as a func- 
tion of the parameters n and Af in the scalar potential 
consistent with Aa/a suggested by the QSO data. A 
few comments are in order at this point. Note that at 
fixed n, ( falls off with increasing M. This behavior fol- 
lows from Eq. (|2.26l) : as the mass scale M increases so 
does (/)o and consequently K,{(f> — 0o)- The dashed line on 
this plot represents the upper bound on the mass scale M 
coming from our initial assumption that the field is track- 
ing today. For a given n, the value of this upper bound 
corresponds to the case when the field has just begun to 
dominate. This is the case of quintessence. Moreover, 
a value larger than this means that the scalar field has 
been dominating the evolution of the universe for a long 
time preventing any formation of structure. Thus every- 
thing above the dashed line is excluded. The equivalence 
principle bound, |^| < 10~^, and the upper bound on M 
produce an allowed band for the pair (n, M) limited by 
the contours —5 and —3 in Fig. ^ 
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FIG. 1: Contour plot for log( with respect to the possible 
choices of n and M in the inverse power law potential. The 
dashed line represents the maximum value for M that allows 
tracking evolution of the field at present. 



The case of quintessence (i.e. along the dashed line in 
Fig. nj is more involved and requires a special approach 
as the field today is no longer in an attractor regime but 
in the transition between the tracker and a scalar field 
dominated solution. Nevertheless, from Fig. ^we expect 
the value of the coupling to be close to C « 10~^. In 
Fig. 121 we show the variation of parameter ^ with n for 
a Universe with fi^ — 0.7 today. This curve was ob- 
tained by solving numerically the evolution equations of 
the field followed by extraction of the coupling C through 
Eq. H2.19|l . One should point out that models with n > 3 
lead to an equation of state for the quintessence field 
in clear disagreement with observations as the latter in- 
dicate 70 = p^/pt/, > 0.4 [s^l- On the other hand, 
models with n < 1 do not alleviate the initial condi- 
tions problem, which was the great advantage of tracking 
quintessence models over the bare cosmological constant. 
Thus, the range of interesting n, 1 < n < 3, in the case of 
quintessence defines the window of the coupling constant 
C that agrees with QSO evidence for changing alpha: 



8.6 X 10"^ < C < 1.4 X 10" 



(2.28) 



We can calculate now the value of d/a at present for 
a given n using Eq. (|2.20|l . In Fig. |3|we show the shapes 
of these curves for the tracker and quintessence cases. 

Finally, as in the case of the exponential potential, the 
BBN bounds on Aa are satisfied with this model. The 
evolution of the field in the radiation dominated epoch is 



«(</) - 0o) = ncbo [-1 + (1 + Zeq)-^/(2+n) 

-4/(2+n)' 



1 + Z 
1 + Zeq 



(2.29) 



FIG. 2: Variation of ^ with the parameter n for an inverse 
power law quintessence model. 



,x 10 




FIG. 3: a/ a at present for a tracker evolution along an inverse 
power law potential (solid line) and for a quintessence model 
(dot-dashed line). 



10"^ < \Aa/a{z = lO^^)] < 4 x 10-^ well within the 
bound of Eq. H1.5|l . 



C. SUGRA type potentials 

Let us now consider the more general potential studied 
in Ref. 1381 which we write in the form 



hence, we can verify that 



(2.30) 



See Refs. |33,|39j that motivate this form within a super- 
gravity context, hence the name adopted here. For nega- 
tive n this potential has a minimum with non-vanishing 
energy density which upon convenient choice of the en- 
ergy scale Vq can explain the present acceleration of the 
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universe. The evolution of the field cannot be written 
in an explicit form but, for [ncj))^ ^ n/\f3 the evolu- 
tion in the tracker regime can be written in terms of a 
perturbative solution to order k written recursively, such 
that 



log ICI 



+ 37ln(l + z) 



ln(K(5 



)ffc-i) 



-2 In 



where 



= In 



2 



7(2-7) 



(2.31) 



(2.32) 



and the last logarithm term only exists if A 7^ 0. The 
case {K(t))^ ^ n/XP is simply the power law case that we 
studied in part B. 

Following the method adopted above to estimate C, its 
value is given by 



10- 



K(l)(k){z = 3) - = 0) ' 



(2.33) 



where (/((j.) is given by Eq. (|2.31|l . In Fig. 0] we show a 
contour plot for the value of C, necessary to account for the 
QSO data, using the second order approximation k = 2. 
In the tracker regime we expect a weak dependence of any 
derived quantities on the parameter n. This argument 
follows from the initial assumption that the field is rolling 
along the exponential side of the potential. We have fixed 
f3 — 2, n — —11 and made A and Vq free parameters. The 
particular case A = 1/2 leads to the form of the potential 
in Ref. [s^ • In Fig.^we show a contour plot for log C with 
respect to the possible choices of n and Vq. The dashed 
line represents the value of Vq that corresponds to the 
present energy density of dark energy assuming that the 
field has already reached the minimum of the potential 
at z = 0. This method of estimating Vq is particulary 
good for A > 1. 

As we have seen in the inverse power law potential case, 
as we approach the dashed line (or quintessence like be- 
havior) the value of the coupling approaches ( ~ 10^^. 
For very small values of Vq and large A we enter the region 
of parameters disfavored by the equivalence principle ex- 
periments. Again, as in the inverse power law potential, 
we obtain an allowed band, now for the pair (A, Vb), lim- 
ited by the contours —5 and —3. For increasing (3 one 
can verify that this band becomes narrower. 

To calculate the variation of a at the time of nucleosyn- 
thesis we just have to realize that the only modification 
in Eq. (|2.31() is that in the radiation era 



X{K(j)f^^^ = + 31n(l zeq) + 41n 



1 



1 + Zcq 



(2.34) 



In Figs. [SI and we show that both the values oi a/a 
at present and Aa/a at nucleosynthesis for this model 
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FIG. 4: Contour plot for log \ (\ with respect to the possible 
choices of A and Vo for the model of Eq. I|2.3fl^ with n = — 11 
and 13 = 2. 



are within the current experimental and observational 
bounds. A striking property is that these quantities are 
extremely insensitive to the values chosen for the param- 
eters of the model, a feature we came across before in 
the previous models. This should not come as a surprise 
because when we estimate C by fixing Aa/a = 10~^ at 
z = 3, we are effectively demanding that the late time 
evolution is roughly the same for any choice of parame- 
ters. 
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FIG. 5: Contour plot for a/ a with respect to the possible 
choices of A and Vo for the model of Eq. 12.301 with n = — 11 
and (3 = 2. 
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FIG. 6: Contour plot for Aa/a at nucleosynthesis with re- 
spect to the possible choices of A and Vo for the model of 
Eq. llOnt with n = -11 and /3 = 2. 



D. Coupled quintessence 

We now consider a linear coupling between the scalar 
field and the background fluid such that Eqns. H2.16I) and 
(|2.17l) are now written in the form: 



Pb + 3Hjpb 



— = kC[pb 



3pb) ■ 



(2.35) 
(2.36) 



Note that the background energy density and pressure 
enter in the combination ps 3pB to account for the 
fact that the background radiation does not provide any 
source for cj). 

We assume an exponential potential for the scalar field 
V{(j)) = Vo exp(— Ak(/)). Such a coupling can arise after a 
conformal transformation to the Einstein frame in Brans- 
Dicke theory 0|. We note, however, that our model 
cannot be considered as pure Brans-Dicke theory as the 
latter does not provide a <j>F^^F^''' coupling. 

It is known that the system has two cosmologically 
interesting attractor solutions 0) \^ '■ 

(1) a scaling solution with 



70 



C'^-CX + 37 
(C- A)2 

C^-CX + 37 



(2.37) 
(2.38) 



(2) a scalar field dominated solution with il^ = 1 and 
70 = /3, identical to the minimally coupled case (see 
Sec. EH). 

For solution (1), the power law expansion a cx i^* is 
given by 



P ■ 



37 V A 



(2.39) 



(not to be confused with pressure). The solution is in- 
fiationary for p > 1, that is, for —A < 2C for a matter 
background. 

The evolution of the field for this kind of coupling is 



37 
X-C 



ln(l -f z) , 



(2.40) 



which reduces to Eq. (|2.21() for C = 0. Therefore, fol- 
lowing the same approach as before, we have that in this 
case the coupling of the scalar field to the electromagnetic 
Lagrangian is 



c 



10- 



A-C 



ln(4) 37 



(2.41) 



in order to explain the Webb et al. variation of a. 

Let us consider the scenario where the scalar field is 
coupled primarily to the dark matter and its coupling to 
baryons is mediated only by photons via c/jF^j^^F^'^ inter- 
action to suppress the violation of equivalence principle. 
In Fig. [7| we show the late time evolution of the energy 
densities of the scalar field, dark matter, baryons and 
radiation in the case when the scalar field is coupled to 
all of the dark matter. We see that at early times the 




2 1 
log (1+z) 



FIG. 7: The late time evolution of the energy densities of the 
scalar field p^, dark matter PdM) baryons pj^i and radiation 
p^, in the coupled quintessence model when all of the dark 
matter is coupled to the scalar field. 

field scales like radiation (as Prad — /Orad/3 and hence the 
right hand side of Eqs. H2.35|l and H2.36|l vanishes), and 
it scales like ordinary matter after the radiation-matter 
transition. At this stage one can take C to be effectively 
zero in Eqs. H2.35|l and (|2.36() following the assumption 
that the (j) couples to baryons very weakly. Only when 
the dark matter contribution becomes important, does 
the scalar field reach the attractor solution (1) above, 
driving the universe to an accelerated expansion with 7^ 
very close to zero, in good agreement with observations. 
However, the dark matter should play an important role 
in the process of structure formation, which would not 
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be possible if it becomes dominant at such a late stage. 
For this scenario to be viable we are forced to introduce 
a component of dark matter that is coupled to the scalar 
field and a component that is not. An alternative ap- 
proach, Ref. jyo\, would be to postulate a non- linear cou- 
pling C = C(<^y. 

In Fig.|Slwe show the evolution of Aa/a for two differ- 
ent values of the slope of the potential A. The damped os- 
cillations in the direction of small redshift reflect the fact 
that the solution is approaching the attractor. In Fig. |^ 
we show the dependence of the coupling ( on the slope of 
the potential A assuming different present contributions 
of the dark matter component coupled to the scalar field. 
The value of C can be extracted from Eq. (|2.37() where 
now we also have to take into account bi-component dark 
matter. Therefore, 

^^^_Wa^£T2^^ (2.42) 

is a good estimate for the value of the coupling in the 
tracker regime. Here flcUM = ^cBu/ i^cUM + ^4>) and 
r^cDM stands for the density parameter of the dark matter 
component coupled to the scalar field, 

We assume that ft^ « 0.7 today. Realistic models of 
coupled quintessence must have a small flcUM for suc- 
cessful structure formation. Consequently, we see from 
Fig.ini that the equivalence principle constraints can only 
be satisfied for larger values of A as the amount of (/)- 
independent component of dark matter increases. For 
large values of A we have seen from Fig. [S] that the at- 
tractor is reached at a later stage, hence the field and 
necessarily Aa/a oscillate heavily at low redshifts. 




z 



FIG. 8: Evolution of Aa/a for a coupled quintessence model 
against the observational QSO data. The coupled dark matter 
component has a contribution JIcDM = 0.05 and Jl^ = 0.7 
today. The solid curve corresponds to A = —100 and the 
dashed curve to A = —10. 

As in the case of a pure exponential potential it turns 
out that a/a = 4.8 x 10~^^ yr~^- However note that this 
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FIG. 9: Dependence of on the parameter A in the cou- 
pled quintessence model. The curves correspond to different 
contributions of the coupled dark. From top to bottom we 
have QcDM = 0.02, 0.05, 0.25, and we assume = 0.7 to- 
day. The dashed curve corresponds to the variation of ^ in 
the uncoupled case of Sec. 12 Al 

is the value that corresponds to the attractor solution. 
Since the field is oscillating around the latter, the value 
of d/a can actually be larger by an order of magnitude 
or smaller, depending on the phase of oscillations. 

The calculation of a at the time of primordial nucle- 
osynthesis is more involved as it requires the calculation 
of the redshift where the attractor behaviour is reached. 
Our estimate for the relative change of a at BBN time, 
Aa/a « 10""^, is consistent with limit (|1.5|l . 

E. Other quintessence models 

In this subsection we analyze some other quintessence 
models suggested in the literature and evaluate if a cou- 
pling to the electromagnetic field can provide a realistic 
scenario compatible with terrestrial and cosmological ob- 
servations. In Fig. ^1 we plot the late time evolution of 
Aa/a for three models, the two exponential case [i^ . 
the Albrecht-Skordis model l48l|, and variations to the 
SUGRA-inspired model proposed in 38]. Also included 
in the plot are the data points of Webb et al . In Ta- 
ble we collect other relevant pieces of information such 
as typical values of the quantities C: a/ct and Aa/a at 
nucleosynthesis. In addition, we indicate which models 
exhibit oscillatory behaviour at late times. 

The SUGRA model (h) in the table is the original one 
of Ref. , in which the field rolls along the power law 
slope of the potential whereas in the remaining cases the 
field starts rolling from the exponential side. This ta- 
ble suggests that various models can be consistent with 
the present day constraints on a/ a and with the BBN 
bounds, however only a few satisfy the Oklo bound and 
only at the cost of severe fine tuning. The latter models 
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Model 


A 


B 


C X 10^ 


a/a X lO^*' 


Aa/a{z = 0.14) 


Aa/a{z = 10'°) 


osc 


2EXP a 


10 


-8 


-20 


8.6 


-3.2 X 10~® 


1.6 X 10"^ 


V 


b 


15 


0.1 


3.9 


4.0 


-8 X 10"^ 


-2.2 X lO""* 


X 


c 


10 


-0.5 


3.9 


-0.2 


-7 X 10"* 


-3.2 X 10"" 


V 


A-S d 


10 




-50 


-1.7 


2 X 10"^ 


3.5 X 10-=* 


V 


e 


6 


lA/A^ 


4.5 


1.2 


-2 X lO"'' 


-5.5 X 10"'' 


X 


f 


6 


0.985/^2 


11.2 


-1.4 


-1 X 10"* 


-1.4 X 10"^ 


V 


g 


8.5 


0.93/yl^ 


-30 


-4.7 


-1 X 10"* 


2.5 X 10"^ 


V 


SUGRA h 


0.5 


-11 


1.08 


4.4 


-9 X lO"'^ 


-0.3 X lO""* 


X 


i 


0.5 


-11 


-0.85 


4.5 


-9 X lO"'^ 


-0.9 X lO""* 


X 


j 


20 


-2 


25 


10.7 


-3 X 10"" 


4.4 X 10"'' 


V 


k 


2.2 


-2 


-1.7 


-0.8 


-2 X 10~* 


-9 X 10"^ 


V 



TABLE I: Approximate values for a/a{z = 0) and Aa/a for BBN (z — 10^") and Oklo phenomenon {z = 0.14) epochs 
for several quintessence models. 2EXP: V = Vb {e^'"*' + e^'^'^) , Ref. lH; A-S: V = K-^e"'*'"* [{K.(j> - Cf + B)] , Ref. |3; 
SUGRA: V = Vbexp {A{K,(l>f) {nfj))'^ , Ref. |3il. In the Sugra model (h) the initial condition of the field is (re^in)^ < B/2A 
and in model (i) {K(j)in)^ 3> B/2A. The models which have late time oscillations of the field have a tick in the last column. We 
have assumed — 0.7 at present. 



are labeled with (c), (f), (g) and (k) in Tabled In these 
models, the field performs a low amplitude late time os- 
cillation hence satisfying the Oklo bound. Models (f ) and 
(g) were recently analysed in Ref. |5l| . The possibility for 
the kind of behaviour of model (g) was pointed out before 
in Ref. 15j. It is clear from Fig. IIUI and Table that in 
order to have a realistic way of constraining quintessence 
models we need to employ the wealth of data available 
on the variations of a over time scales much larger than 
those observed through the QSO observations. 




z 



FIG. 10: Late time evolution of Aa/a for the models in Table 
13 that pass the Oklo bound: 2EXP model (c) in dashed line; 
A-S models (f) and (g) of Ref. p3l in light solid and solid lines 
respectively; SUGRA model (k) in dot-dashed line. 

In the next section we propose a mechanism that en- 
able us to reconcile a large detection of a variation on a 



at large redshift with the tight Oklo constraint without 
need of fine tuning in the parameters of the potential. 

3. NON-MINIMAL COUPLING OF 
QUINTESSENCE TO THE ELECTROMAGNETIC 
FIELD 

In this section, we propose a radical alternative to the 
allowed couplings between the quintessence and electro- 
magnetic field, which will allow us to reconcile all the 
current constraints on the variation of a with the data. 
We start by asking what motivates a linear coupling be- 
tween quintessence and the electromagnetic field? Obvi- 
ously, there are two simple arguments for this choice: its 
apparent simplicity and the fact that almost any inter- 
action can be reduced to this form if Bp{(l)) can be ex- 
panded as a series in (j){z) — 0o provided that the latter is 
small. This linearization can break down for larger (/'(z), 
and thus the comparison with Aa/a at very high z (i.e. 
BBN limits) can only tell us whether the linear depen- 
dence needs to be modified. Of course the linear choice is 
not the only possibility, an d quadr atic couplings between 
matter and the scalar field 13, 50] and const-fexp {—Kcf)} 
couplings [3| being also actively discussed in the past. 
All these approaches share the same weakness: they do 
not calculate these dependences from first principles but 
rather simply postulate them. In many ways it is anal- 
ogous to simply choosing a potential for a tracker field 
by hand. Clearly, when an ad hoc choice of Bp[(l)) is al- 
lowed, the profile of a{t) cannot be calculated and any 
a{t) is possible. In particular, one can "engineer" such a 
function Bp{(j)) such that both the Oklo constraint and 
Webb et al. data points are satisfied. Moreover, one 
can choose dBp{(j)) / dMz = 0) <C k, suppressing the fifth 
force mediated by cj) |2Jj HOl ■ 
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We would like to point out another rather generic sce- 
nario where the effective fine structure constant is mod- 
ified in atomic physics but has almost no impact on nu- 
clear and particle physics. This would allow us to rec- 
oncile the Oklo constraints with the results of Webb et 
al. and suppress the non-universality of the gravitational 
force mediated by (/>. To do this we propose a momentum 
dependent modification of the interaction Hl.l|l . Quite 
generically, the vertex between two photons and a (j) field 
can be written in the momentum representation as 



1 



1 



+ 9f^i^ P 



I[p^,q^, (pq),' 



(3.43) 



where p ± q/2 arc photon momenta and q is the mo- 
mentum of the scalar field. Eq. (|3.43|) has introduced 
a form factor / that depends on the momenta and the 
value of the scalar field. The minimal choice of coupling, 
/[p^, , (pq), 0] = i^Cj discussed in previous sections. 
For all practical purposes, one can neglect the momen- 
tum of the scalar field so that 



f[p',qMpq)] = f[p^ 



(3.44) 



Let us now imagine that this coupling has a form factor 
which tends to zero as becomes larger than a certain 
critical momentum scale fi. For definiteness, we assume 
a power-like suppression, 



f[p'] = «c- 



(3.45) 



where ^ is the momentum independent fundamental con- 
stant. This choice of / leads to the anomalous depen- 
dence of the effective coupling constant on the photon 
virtuality such that 



Ccff(/) - c 



a{p^) = ac{l + KCcff0) with 
_ / C for < ^2 

^2 _p2 I Q fgj. |p2| ^ ^2 



(3.46) 



In other words, the high virtuality photons always couple 
to matter with a constant time-independent strength ac, 
while the couplings of low virtuality photons are modi- 
fied by the presence of 4> which evolves with cosmological 
time. Of course, ac will have a usual high-energy log- 
type running with the momentum. In the language of 
operators the choice ()3.45|) would correspond to a new 
non-local operator that can be written as 



1 



-4>F 



Q2 



(3.47) 



The zeroth order term in the expansion over j \j? cor- 
responds to a usual linearized form of H2.13|l that we 
were using in the previous sections. We also note that 



Eq. (|3.47|) does not allow for arbitrary additive redifi- 
nition of at z = 0. As before, we assume that the 
cosmological evolution of is driven by its self potential. 

In the absence of the form factor, the effective cou- 
pling between nucleons and scalar field, (pfK(j)NN (see 
the notations of Ref. |0|), induced via the coupling to 
photons is given by Cat ~ 10~^C < 10~^ [sj, leading to 
the bound (|1.7|l . It is clear that in order to obtain a fur- 
ther suppression of coupling by the form factor (|3.45() 
one should choose the scale p, to be lower than the char- 
acteristic virtualities of photons involved in the hadronic 
and nuclear matrix elements. The latter are of the order 
of the inverse characteristic size of a system, which in 
the case of coupling is the size of the nucleon. For 
example, the coupling of (f> to nucleons. Cat, will receive a 
typical suppression of ^2/(200 MeV)^ for fi < 200 MeV: 



(200 MeV)2 



(3.48) 



where 200 MeV is roughly the inverse size of a nucleon. 
This suppression can be even stronger for a steeper form 
factor. This leads to a considerable weakening of all 
constraints coming from the non-universality due to the 
(/)-exchange. For example, choosing /i around 10 MeV, 
we get an additional 10^^ — 10^^ suppression so that 
the effective coupling to nucleons is now much smaller 
Cn ^ (10~^ — 10~^)C. This results in a much weaker, 
C ^ 1, bound than eq. ((r7|l . 

Likewise, the bounds coming from Oklo phenomenon 
will also be weakened. This is because in the presence 
of the form factor (|3.46|l the dependence of nuclear lev- 
els on (/) entering in the standard calculation 0, 113 
is suppressed and consequently these levels are less af- 
fected by the change in the (j) field. Indeed, the ex- 
pected weakening of the Oklo bound is at least a factor 
of Ai^/(50 MeV)^ ~ 0.04 with the same choice of ^. Here 
50 MeV is a typical inverse radius of a heavy nucleus. 
This is just enough to make the results of Oklo analysis 
H m iHOl compatible with Webb et al. even with d ~ 
const and without any fine tuning. 

On the other hand, the atomic physics phenomena 
for which the characteristic photon virtualities are much 
smaller than ^ ~ 10 MeV will be sensitive to "dressed" a 
that evolves as a function of (p. In particular, this refers 
to the relativistic corrections to the atomic energy lev- 
els, used for the QSO analysis of Webb et al. Thus, the 
existence of the form factor allows to "isolate" the phe- 
nomenon of changing alpha to atomic physics and remove 
it from nuclear and particle physics. 

The anomalous dependence of the fine structure con- 
stant on the photon momentum (|3.46() can be probed 
with high precision QED measurements, as different ex- 
periments probe different photon virtualities. Some of 
these measurements are performed with atomic and con- 
densed matter systems, where the characteristic momen- 
tum transfers are much smaller than fj,. Take, for exam- 
ple, the determination of a from the quantum Hall effect 
where the characteristic momenta of photons are very 
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small. This experiment will be sensitive to the "dressed" 
value of a(jP ~ 0) with obvious condition (p = (j){z = 0). 
On the other hand, the precision measurements sensitive 
to QED radiative corrections such as electron and muon 
anomalous magnetic moments can also probe a at vir- 
tualites comparable to fi where the value of a is given 
by H3.45|l . The expression for the magnetic dipole mo- 
ment of the electron and muon including the first loop 
correction will be given by 



e(p2 ~ 0) 



nii 



1 



a{p^ ~ mf) 
2^ 



(3.49) 



where i — e,^. The electron g — 2 experiment |52| has 
measured the anomalous magnetic moment of the elec- 
tron to an accuracy that allows for the extraction of 1 /a 
to nine digits. Now we have to take into account that the 
one-loop correction has some photons with virtualities 
^ ^"^ ^ (10 MeV)^. This will introduce a correction 
to 5 — 2 of the electron due to the anomalous behaviour 
of a at large momenta at the level of 



A(5-2), 



\2ttJ [ u? 



10"^"- 10 



-11 



(3.50) 



for ~ 10 MeV and typical C,K(f) ~ 10^'"' suggested by 
Webb et al. results. This correction is taken relative to 
the value of g — 2 calculated with "fully dressed" value of 
a (|3.46|l . Since (g— 2)e is used for the extraction of a, one 
should absorb correction (|3.50|l into the measured value 
of a and compare it with the second best determination 
of the fine structure constant from the quantum Hall ef- 
fect. The correction (I3.50|l corresponds to the shift of l/a 
in the eighth digit which appears to be marginally con- 
sistent with the determination of a from the solid state 
physics 01 • 

The largest effect is expected in the g — 2 of the muon, 
where the radiative corrections probe "undressed" a, or 
Q!c, (|3.46() at virtualities comparable to the muon mass, 
which is much larger than a chosen value of /i ~ 10 MeV. 
Therefore, the expected correction to the muon anoma- 
lous magnetic moment (relative to the value one would 
expect with "dressed" a given by will be on the 

order 



A(5-2), 



(3.51) 



which is larger than both the experimental and theoret- 
ical accuracy of [g — 2)^, unless C,K(f) < 10~^, where (j) 
is the value of the scalar field now. This might still be 
marginally consistent with Ref. '5| if (K(j){z ~ 0.5) ~ 
10-^ 

Therefore, the anomalous running of a with the mo- 
mentum H3.46|l . H3.46|l comes tantalizingly close to the 
accuracy of modern tests of QED. At the same time, it 
allows us to disconnect the phenomenon of "changing a" 
from nuclear and hadronic physics. This way the results 



of Webb et al. coming from the atomic physics, are not 
in contradiction with Oklo phenomenon and are of no 
immediate consequence for the fifth- force experiments. 



4. CONCLUSIONS 

We have studied a variety of cosmological scalar field 
models coupled to the electromagnetic Lagrangian. The 
evolution of the scalar field over cosmological times in- 
duces the effective change of the fine structure constant. 
We have considered two generic possibilities: 

(i) Tracker solution. We have shown that for a wide 
range of models where the scalar field is tracking, there 
exists a region in parameter space which allows concor- 
dance with the QSO data on non-zero Aa/a ^ with 
other constraints on Aa/a coming from the Big Bang 
nucleosynthesis as well as the precision measurements of 
the universality of the gravitational force. The predic- 
tions for the present day a/ a are typically within one 
order of magnitude from current experiments and within 
the reach of near future laboratory experiments. Typical 
constraints on the parameter space are shown in Figs. 
and^ There is an upper bound on the mass scale M or 
Vq from the requirement that the tracking solution per- 
sists until z = and <ti Pg today. The consequence 
of this is the lower bound on C, the coupling between 
the scalar field and the electromagnetic F^^F^^'^ . This 
typically leads to Q > 10^^ (see Figs. and from the 
dashed line in Fig.lSJ in the region of parameters studied. 
The universality of the gravitational force puts an upper 
bound on the couphng, C, < 10"'^. 

(ii) Dark energy. We have found that the simple in- 
verse power law model for quintessence also offers a small 
region of parameter space in agreement with the QSO 
observations. Here, for a chosen n the mass scale M is 
fixed, therefore, the only free parameter remaining is the 
coupling To fit the QSO data for a chosen n we obtain 
one single Q (or a small window due to the uncertain- 
ties related with the data). We have also considered the 
model of the scalar field coupled to the dark matter. We 
have found that only large values (if A is negative) of the 
slope of the scalar potential allow this model to pass the 
equivalence principle constraints. This model can gen- 
erate oscillatory behaviour for a{z) which may lead to 
the enhancement of the present day signal of a /a. We 
looked at other models of quintessence suggested in the 
literature. Their relevant properties for this work are 
summarised in Table ^ We expect that complimentary 
information on the present equation of state of the uni- 
verse, on the coupling C, and on the ratio a/ a can help 
us to back up or rule out models of quintessence and sce- 
narios in which a dynamical subdominant scalar field is 
evolving in a tracking regime. 

A notable feature of this detailed analysis is that 
apart from a small region of parameter space in some 
quintessence models with late time oscillations (e.g., 
models (c), (f), (g) and (k) of Table P) we have not 
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come accross a model that satisfies both the Oklo bound 

|Aa/a| < 10""^ at z ~ 0.14 and Aa/a 10"^ at z - 1 

reported by Webb et at. In an attempt to reconcile this, 
we have suggested that the existence of form factor in 
the coupling of C at around 10 MeV with respect to the 
photon momentum may lead to the suppression of the 
effective coupling between (f> and nucleons and nuclei. 
This is sufficient to isolate the phenomenon of changing 
couplings to the realm of atomic physics, while relaxing 
the bounds coming from the Oklo phenomenon and the 
fifth- force experiments. This amounts to postulating an 
anomalous running of a with momentum and thus gives 



an additional possibility of checking this idea with high- 
precision QED measurments, sensitive to various ranges 
of the photon virtualities. 
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